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Abstract: By using the complete discrimination system for polynomials, we study 
the number of positive solutions in C[0, 1] to the integral equation 
P( x ) = Jo k ( x ,y)<fi"(y) d y, where k{x,y) = ipi(x)(j>i_(y) + ^(x^iy), (p t (x) > 
0,<f>i(y) > 0, < x,y < l,i = 1,2, are continuous functions on [0,1], n is 
a positive integer. We prove the following results: when n = 1, either 
there does not exist, or there exist infinitely many positive solutions in 
C[0, 1]; when n > 2, there exist at least 1, at most n + 1 positive solutions 
in C[0, 1]. Necessary and sufficient conditions are derived for the cases: 
1) n = 1, there exist positive solutions; 2) n > 2, there exist exactly 
m(m € {1, 2, ■ ■ • , n + 1}) positive solutions. Our results generalize the ex- 
isting results in the literature, and their usefulness is shown by examples 
presented in this paper. 

Keywords: Integral Equations, Positive Solutions, the Complete Discrimination Sys- 
tem for Polynomials, the Number of Solutions 



1 Introduction 



The existence of positive solutions to integral equations is an active research field and has 
important applications in the stability of feedback systems [1,2]. In 1991, the number of positive 
solutions to the following integral equation 
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ip(x) = [ k(x,y)ip 2 (y)dy (1) 
Jo 

was discussed in [3]. In this paper, we will study the number of positive solutions in C[0, 1] to 
following more general integral equation 

<p(x) = [ k{x,y)<p n {y)dy (2) 
Jo 

where 

k{x,y) = tp 1 (x)(j) 1 (y) + <p 2 {x)<t>2{y),ipi{x) > Q,4>i{y) > 0,0 < x,y < l,i = 1,2 

are continuous functions on [0, 1], n is a positive integer. We prove the following results: when 
n = 1, either there does not exist, or there exist infinitely many positive solutions in C[0, 1]; 
when n > 2, there exist at least 1, at most n + 1 positive solutions in C[0, 1]. Especially, when 
n is an odd number greater than 2, there exist at least 1, at most n positive solutions in C[0, 1]. 
Necessary and sufficient conditions are derived for the cases: 1) n = 1, there exist positive 
solutions in C[0, 1]; 2) n > 2, there exist exactly m(m e {1,2, • • • ,n + 1}) positive solutions 
in C[0, 1]. Our results generalize the existing results in the literature, and their usefulness is 
shown by examples presented in this paper. 

In essence, the number of positive solutions to (2) can be transformed into determination 
of real roots of a certain polynomial, which is a century- long, albeit still active research area in 
mathematics. The classical Sturm method or Newton formula can be employed to determine the 
real root distribution of polynomials [4-7], but the Sturm method is inefficient in establishing 
discriminant systems for high-order polynomials with symbolic coefficients [4,6,7], and Newton 
formula involves in recursive procedure to determine the real roots, thus it is difficult to establish 
explicit criteria [4,5,7]. 

More recently, Yang et al. established the complete discrimination system for polynomials, 
which can give a set of explicit expressions based on the coefficients of polynomials to determine 
the root distribution of polynomials [6,7]. 

Let 

f(x) = a x n + a x x n - x + ■ ■ ■ + a n e P n , 

the Sylvester matrix of f(x) and its derivative / (a;) [6,7] 

ao a i a 2 ■ ■ ■ dn-l &n 
na (n - l)ai ■ ■ ■ 2a„_ 2 a„_i 

a a l ' ' ' a n~2 ffln-1 

na ■ ■ ■ 3a„_ 3 2a„_ 2 

a ai a 2 • • • a„ 
na (n — l)ai • • • a„_i 

is called the discrimination matrix of f(x), denoted as Discr(f). 

[£>i(/),£> 2 (/),•• •,£>„(/)] 

the even-order principal minor sequence of Discr(f), is called the discriminant sequence of 
fix). 



a n -i 
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[sign(Di), sign(D 2 ), ■■■ , sign(D n )} 

is called the sign list of the discriminant sequence [Di,D 2 , • ■ ■ , D n ], where sign(-) is the sign 
function, i.e., 



1 


if x > 0, 





if x = 0, 


-1 


if x < 0. 



Given a sign list [s\, s 2 , • • • , s n ], we can construct a revised sign list 

[ei)£2) • • • j £n] 

as follows: 

l)If [s{, Sj+i, • • • , Sj+j] is a section of the given sign list and Sj 7^ 0; Sj+i = Si+2 = ■ • • = 
Si+j-i = 0; Si+j 7^ 0, then replace the subsection consisting of all elements 

[Sj+l, Sj+2, • • • , Sj+j-i] 

by the following subsection with equal number of terms 

[ s i: s i: s i: s i, Si, Si, Si, Si, Si, ' ' '] 



i.e.,£ i+r = (-1) 2 . St ,r = 1,2, ••• ,j - 1. 

2)Let £fe = Sk for all other terms, i.e., all other terms remain the same. 

Lemma 1 [6,7] Given the polynomial with real coefficients f(x) = a^x n +aix n ^ 1 +■ ■ -+a n <E 
P n . If the number of sign changes in the revised sign list of its discriminant sequence is v, and 
the number of non-zero elements in the revised sign list is /1, then the number of distinct real 
roots of f(x) is fi — 2v. 

Remark 1 The discriminant sequence of f(x) can also be constructed by the principal 
minors of the Bczout matrix of f(x) and / (x) [6,7]; the number of distinct real roots of f(x) 
can also be determined by the sign difference of Bezout matrix of /(x) and / (x) [6,7]. 

Remark 2 The complete discrimination system for polynomials can also be used to deter- 
mine the number and the multiplicity of complex roots [6,7]. 

Yang and Xia also proposed a method to determine the number of positive (negative) roots 
of a polynomial [8], which is similar to Lemma 1 in principle, but is more efficient. 

Lemma 2 [8] Given the polynomial with real coefficients f(x) = a^x n +aix"~ 1 +• • - + a n G 
P n , a 7^ 0, a n ^ 0. Let h(x) = /(— x) and {d\,d,2, • ■ ■ , c?2n+i} be the sequence of the principal 
minors of the discrimination matrix Discr(h) of h(x). If the number of sign changes in the 
revised sign list of the sequence {^1(^2,^2^3, • • • , d2nC?2n+i} is v, and the number of non-zero 
elements in the revised sign list is \i, then the number of distinct positive roots of /(x) is fi — 2u. 



2 Main Results 



Consider the problem of determining the number of positive solutions in C[0, 1] to the integral 
equation of the following form 



/' 

Jo 



<p(x) - / k{x,y)<p n {y)dy (3) 
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where 



k{x,y) = ip 1 (x)(j) 1 (y) + ip 2 (x)(j) 2 (y) 7 ip i (x) > 0,<f>i(y) > 0, < x, y < 1, i = 1, 2 

are continuous functions on [0, l],n is a positive integer. 
Denote 

a n - iii = C l n <j)i{y)ipl~\y)ip l 2 {y)dy, i = 0,l,---,n, 
Jo 

b n - i ,i = C i n [ Mv)<Pi~ i (v)<P2(v)dV> i = 0,l,---,n, 
Jo 

Q!» = &n-*,i — ttn-i+ljj-l, i = 1,2, •■•,n, ao = 6 n ,o, Ck n+ i = — ao, n , 

where C* , i = 0, 1, • ■ • , n, stand for the combinatorial number. Our main result is as follows: 
Theorem 1 When n = 1, either there does not exist, or there exist infinitely many positive 

solutions in C[0, 1] to the integral equation (3). The necessary and sufficient conditions for the 

existence of positive solutions in C[0, 1] are 01,0 — 1 < and (ai,o — l)(&o,i — 1) — a o, i&i,o = 0. 
Theorem 2 When n > 2, there exist at least 1, at most n + 1 positive solutions in C[0, 1] 

to the integral equation (3). Especially, when n is an odd number greater than 2, there exist 

at least 1, at most n positive solutions in C[0, 1]. 

Theorem 3 When n > 2, the necessary and sufficient conditions for the existence of 

exactly m(m G {1,2, •••,n+ 1}) positive solutions in C[0, 1] to the integral equation (3) 

arc: the number of sign changes v in the revised sign list of the discriminant sequence of 

n+1 

the polynomial f(s) := ctiS 2 ( n+1 ~ 1 ^ and the number of its non-zero elements fi satisfy 

i=0 

m = ' ; or, equivalently, the number of sign changes v in the revised sign list of the sequence 
\d\di, d<2,dz, ■ ■ ■ , din+idin+i} an d the number of its non-zero elements \x satisfy m = v , 
where {d\,d 2 , • • • , d 2n +3} is the sequence of the principal minors of the discriminant matrix 

n+1 

DiscrQi) oih{s) := £ a l (-s) n+1 -\ 

i=0 

Specifically, when n = 2, denote 

ai a 2 a 3 
P = — ,r = — ,t = — < 0, 
a a a 

Ai = p 2 - 3r, A 2 = rp 2 + 3tp - 4r 2 , 
A 3 = -4r 3 + ISrtp + p 2 r 2 - Ap 3 t - 27t 2 , 
[D U D 2 , D 3 , D 4 , D 5 , D 6 ] = [1, -p, -pAi, AiA 2) A 2 A 3 , -tA§]. (4) 

then we have 

Corollary 1 There exist at least 1, at most 3 positive solutions in C[0, 1] to the integral 
equation (1). 

Corollary 2 The necessary and sufficient conditions for the integral equation (1) to have 
exactly 3 positive solutions in C[0, 1] are p < 0, Ai > 0, A 2 > 0, A3 > 0. 

Corollary 3 The necessary and sufficient conditions for the integral equation (1) to have 
exactly 2 positive solutions in C[0, 1] arc p < 0, Ai > 0, A 2 > 0, A 3 = 0. 

Corollary 4 The necessary and sufficient conditions for the integral equation (1) to have 
exactly 1 positive solutions in C[0, 1] are p > 0, or Ai < 0, or A 2 < 0, or A3 < 0. 

Remark 3 If n is even, the integral equation (3) docs not have any negative solutions in 
C[0,1]. 
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Remark 4 If n is odd, since ip(x) is a positive solution in C[0, 1] to the integral equation 
(3) if and only if — ip(x) is a negative solution in C[0, 1] to the integral equation (3), thus, when 
n=l, the integral equation (3) either does not have, or has infinitely many negative solutions 
in C[0, 1]; when n is odd and greater than 2, the integral equation (3) has at least 1, at most 
n negative solutions in C[0, 1]. 

Remark 5 When n = 1 , the necessary and sufficient conditions for existence of negative 
solutions in C[0, 1] to the integral equation (3) are the same as in Theorem 1. When n is 
odd and greater than 2, the necessary and sufficient conditions for existence of exactly m(m G 
{1, 2, • • • , n}) negative solutions in C[0, 1] to the integral equation (3) are the same as in Theorem 
3. 

Remark 6 Our method can be extended to the case when the integral kernel k(x,y) is 
i 

taken as ^ ifi(x)4>i(y), where <fi(x) > 0, 4>i(y) > 0, < x, y < 1, i = 1, 2, • • • , I, are continuous 
functions on [0, 1]. 

Remark 7 The conclusions in [3] are equivalent to Corollaries 1,2,3 above. 



3 Proof of the Theorems 



Proof of Theorem 1 When n = 1, the integral equation (3) becomes 

ip(x) = [ k(x,y)<fi(y)dy (5) 



Thus, we have 

<p{x) = (pi(x) / (j) 1 {y)ip{y)dy + ip2{x) \ <p 2 {y)(p{y)dy 
Jo Jo 

If ip(x) is a positive solution in C[0, 1] to equation (5), then ip(x) can be expressed as (p(x) = 

Xiipi(x) + \2(fi2(x), where Ai > 0,A2 > arc coefficients to be determined. Taking it into 

equation (5), we get the following system of algebraic equations 

a l,(Al + a 0,l-^2 = Ai 

hflXi + &o,iA 2 = A 2 

where a lfi = J 4>i(y)<pi(y)dy, o ,i = J <i>i{y)V2(y)dy, b lfi = f fc{y)<Pi{y)dy, & ,i = 
J <^2(y) ( fi2(y)dy. Apparently, the necessary and sufficient conditions for the system of algebraic 
equations (6) to have positive solutions Ai, A 2 are ai.o — 1 < 0, and (ai,o — l)(^o,i — 1) — fl o,i^i,o = 
0. Moreover, if ip(x) is a positive solution to equation (5), then, obviously, for any positive 
constant number c,cip(x) is also a positive solution to equation (5). Thus, there are infinitely 
many positive solutions in C[0, 1] to equation (5). This completes the proof. 
Lemma 3 The system of equations 

a n .ox n + a„_i.ia;™" 1 j/ + a n - 2 .2X n ~ 2 y 2 H h ax.n^xy^ 1 + a , n y n = x , . 

b nS) x n + b n ^ x x n - x y + 6„- 2 , 2 z"-V + ' • • + b x . n ^xy n - x + b , n y n =y ( ' > 

a n -i,i>0, b n - it i>0, i = 0, 1,2, • • • ,n. 
has at least 1, at most n + 1 (at most n, when n is odd) positive solutions, where n > 2. 
Proof Let 

P{x,y) = — — — — -^r— , i — , Q(x,v)- 



a n ,ox n + a n -i^x n 1 y-\ h a ,„y" b nfi x n + b n -i^x n 1 y H 1- b 0>n y r ' 
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x > 0, y > 0, 

then 

1 1 

p(nx, ny) = —— [ p{x, y), q(nx, ny) = —— [ q{x, y), k > 0. 

Let 

E = {x\p(x,l) = q(x,l)}, 

then the number of positive solutions to the system of equations (7) is equal to the number of 
elements in E. In fact, if (x,y) is a positive solution to (7), then 

P{x,y) = q(x,y) = 1, 

thus | e E. Conversely, if x 6 E, since n > 2, it is easy to verify that ( n ~y/p(x, "~-^/p(x, 1)) 

is a positive solution to (7). 

Suppose x E E, by p(a;, 1) = q(x, 1), we have 

b n flx n+1 + (6 n _i,i - a n ,o)a;™ + (6 n -2,2 - a„_i,i)a; n_1 H h (6o, n - ai,n-i)z - ao,n = 

Namely 

ao^™ +1 + Q.\x n + ct2X n ~ 1 + • • • + a n x + a„+i = (8) 

where 

Q-i = b n -i t i - a n _i+i,i_i, i = 1, 2, • • • , n, a = 6„,o, a»+i = -a ,„, 

Since ao = &rt,o > 0, a„+i = — ao,n < 0, equation (8) has at least 1, at most n + 1 positive 
roots. Especially, when n > 2 and is odd, since equation (8) has at least 1 negative root, it has 
at most n positive roots. This completes the proof. 
Proof of Theorems 2 and 3 When n > 2, since 

<fi( x ) = Jo H x iv)v n {y)dy 

= <pi( x ) Jo My)v n (y)dy + ^(x) J^ 1 My) l f n (y)dy 

similar to the proof of Theorem 1, the positive solution ip(x) in C[0, 1] to the integral equation 
(3) can be expressed as <p(x) = Ai<^i(a;) + \2<fi2(x), where Ai > 0, A 2 > are coefficients to be 
determined. Taking it into equation (3), by a simple but lengthy calculation, we see that Ai, A2 
should be positive solutions to the following system of algebraic equations 



(9) 



where 



a n,oA™ + fln-i.iA-L A2 + a„_2,2A 1 A 2 + • ■ ■ + ai !n _iAiA 2 + ao t n 

Xi = Ax 

6„,oA? + 6„-i ; iAr 1 A 2 + &„-2, 2 Ar 2 A| + • • • + fei^-iAiA^ -1 + b , n \% = A 2 

a n -i,i = c i n [ My^i^ivMiy^y, i = o,i,---,n, 

Jo 

bn^^Clf <h(v)<Pi~ i (v)<Pttv)dV, i = 0,l,---,n. 
Jo 

By Lemma 3, we complete the proof of Theorem 2. 

Moreover, from the proof of Lemma 3, we know that finding the positive solutions to the 
system of algebraic equations (9) or (7) can be transformed into finding the positive solutions 
to equation (8). Applying Lemmas 1 and 2 to equation (8), we complete the proof of Theorem 
3. 

Proof of Corollaries 1,2,3,4 Some notations in this proof are defined in Section 2. 
Corollary 1 is a direct consequence of Theorem 2. 
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When n = 2, equation (8) becomes 

a x 3 + aix 2 + a 2 x + a 3 = (10) 

By a direct computation, we know that the discriminant sequence [Di, D 2l D 3 , D4, D5, Dq] of 
the polynomial f(s) := a s 6 + ais 4 + a 2 s 2 + a 3 is determined by (4) (up to a positive factor). 
Since t < 0, it is easy to see that the number of sign changes v in the revised sign list of 
[Di, D 2 , D 3 , D4, D 5 , D 6 ] and the number of its non-zero elements \i satisfy 6 = yU — 2v if and 
only if the revised sign list of [D\, D 2 , D 3 , D4, D 5 , D 6 ] is [1, 1, 1, 1, 1, 1], which is equivalent to 
p < 0, Ai > 0, A 2 > 0, A 3 > 0. This completes the proof of Corollary 2. 

Similarly, the number of sign changes v in the revised sign list of [D\, D 2 , D 3 , D4, D 5 , D 6 ] and 
the number of its non-zero elements /i satisfy 4 = fi — 2v if and only if the revised sign list of 
[Dx,D 2 , D 3 ,D 4 , D 5 ,D 6 ] is [1, 1, 1, 1,0,0], which is equivalent to p < 0, Ai > 0, A 2 > 0, A 3 = 0. 
This completes the proof of Corollary 3. 
Combining Corollaries 1,2,3, we get Corollary 4. 



4 Some Illustrative Examples 



Example 1 Consider the integral equation 

f 6 3 

fix) = {-xy+ -y)<p(y)dy, 0<x<l (11) 

Let 

6 3 

<Pi(x) = gX, <t>i(y) = y, (p 2 {x) = -, 4> 2 {y) = y. 

Then, it is easy to get a 10 = 61,0 = §1 a o,i = &o,i = §■ The conditions in Theorem 1 are met. 
Hence, there are infinitely many positive solutions in C[0, 1]. In fact, ip(x) — c(|x + |), Vc > 
are such solutions. 

Remark 8 From the proof of theorems and the example above, we can see that, for a given 
integral equation, we can not only determine the number of its positive solutions, but also find 
the positive solutions explicitly by solving the algebraic equation (8). 

Example 2 Consider the integral equation 



ip(x) = Jq[18 max{e, — 2x + 1 + e} + max{|, ^(2x — 1 + e)} 
xmax{6,272y-l30}]ip n {y)dy, < x < 1, 

where e > 0, n = 1 or 2. 
Let 

ipi(x) = max{e, — 2x + 1 + e}, <fii(y) = 18, 
<p 2 (x) = max{|, ^{2x - 1 + e)}, 4> 2 {y) = max{6, 272y - 130}. 
When n = 1, using the notations in Section 2 and by a simple computation, we can get 

9 r 1 3 



(12) 



«i,o = / <f>i(y)<Pi(y)dy = 18e+ -, a ,i = J <j>i{y)<p 2 {y)dy = 6e + -; 



h,o = J <f> 2 (y)<Pi(y)dy = 40e+ |, b ,i = J <j) 2 {y)(p 2 {y)dy = 



40 145 
T e + IF 
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Since ai.o — 1 > 0, the conditions in Theorem 1 are not met. Thus, equation (12) does not have 
any positive solutions in C[0, 1]. 

Similarly, when n — 2, by a simple computation, we can get 

a 2 ,o = 3 + 9e + 18e 2 , a M = 6e + 12e 2 , a , 2 = \+ £ + 2e2 ^ 

, n , n 2 , 154 80 2 , „ 145 40 2 
b 2 ,o = 1 + 3e + 40e 2 , o M = —e + —e 2 , b , 2 = 2 + — e + —e 2 ; 

26 73 

ao = h,a = 1 + 3e + 40e 2 , a x = &i,i - a 2j0 = -z-£ 2 + - 3, 

o y 

68 2 17 „ 2 1 

a 2 = o ,2 - ai,i = — — e - — e + 2, 0:3 = -a 0:2 = -2e - e - - . 



Let 



P = — > r= — , t= — ; 

a a-o a-o 

„2 o„ A „2 1 Or ,1^2 



Ai = v - 3r, A 2 = rjr + 3tp - 4r' i 
A 3 = -4r 3 + 18rtp + p 2 r 2 - 4pH - 27t 2 , 
we have (up to a positive factor) 

p= -l + 2.7037e + 2.8889e 2 , 

Ai = -65.959e 2 + 94.716e 3 - 21.593e + 327.26e 4 + 1, 
A 2 = 702. 22e 4 - 1291. 9e 5 + 255. 78e 3 - 2356. 3e 6 - 78.517e 2 - 26.207e + 1, 
A 3 = 1 - 27.778e - 1.4371 x 10 5 e 6 - 23275.0e 4 - 1.0374 x 10 5 e 5 - 63.724e 2 - 1222.6e 3 . 
By numerical computations, it is easy to get 

The real roots of p = arc -1.2197 and 0.2838; 
The real roots of A 1 = are 0.041426 and 0.45024; 
The real roots of A 2 = are -0.70495 and 0.034952; 
The real roots of A 3 = arc -0.21287 and 0.03143. 

Hence, by Corollaries 2,3,4, it is easy to know that there exists a positive number r w 0.03143 
(note here the difference between the exactness of the conditions in Corollaries 2,3,4 and the 
inexactness of the numerical computations above), such that: when < s < ro, equation (12) 
has 3 positive solutions in C[0, 1]; when e = ro, equation (12) has 2 positive solutions in C[0, 1]; 
when e > ro, equation (12) has 1 positive solutions in C[0, 1]. 

Remark 9 The case when n = 2 in the example above has also been studied in [3]. Our 
result here is completely consistent with the result in [3]. 

Example 3 Consider the integral equation (12) in the example above. When n = 3, e is 2 
or 0.2, determine the number of its positive solutions in C[0, 1]. 

Similar to Example 2, when n = 3, using the notations in Section 2 and by a simple computation, 
we can get 

9 27 27 

03,0 = - + 9e + —s 2 + 18e 3 , a 2 ,i = 3e + —s 2 + 18e 3 , 

ai,2 = tj£ 2 + 6e 3 + e, a ,3 = |e 3 + + ^ e + 

3 9 163 

03,0 = - + 3e + -e 2 + 40e 3 , 62,1 = e + — e 2 + 40e 3 , 
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7 + 3e + ^e 2 + 40e 3 , «i = -8e + ^-e 2 + 22e 3 - -, 
4 2 3 4 




Hence, when e = 2, we have a s s + a lS 6 + a 2 s 4 + a 3 s 2 + a 4 = 344.75s 8 + 212.42s 6 - 18.889s 4 - 
40.431s 2 — 8.0833. By a simple computation, the revised sign list of its discriminant sequence is 



By Theorem 3, equation (12) has only 1 positive solution in C[0, 1]. 

When e = 0.2, we have a s 8 + a lS 6 + a 2 s i + a 3 s 2 + a 4 = 1.85s 8 - 3.1273s 6 + 0.68711s 4 + 
0.63519s 2 — 0.17533. By a simple computation, the revised sign list of its discriminant sequence 
is 
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[1,-1,-1,-1,1,1,1,-1]. 
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